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The anomalous scaling behavior of the topological susceptibility Xt m two-dimensional 
CP 1 sigma models for N < 3 is studied using the overlap Dirac operator construction of 
the lattice topological charge density. The divergence of \t hi these models is traced to the 
presence of small instantons with a radius of order a (= lattice spacing), which are directly 
observed on the lattice. The observation of these small instantons provides detailed confirma- 
tion of Liischer's argument that such short-distance excitations, with quantized topological 
charge, should be the dominant topological fluctuations in CP 1 and CP 2 , leading to a diver- 
gent topological susceptibility in the continuum limit. For the CP N ~ X models with N > 3 
the topological susceptibility is observed to scale properly with the mass gap. These larger 
N models are not dominated by instantons, but rather by coherent, one-dimensional regions 
of topological charge which can be interpreted as domain wall or Wilson line excitations and 
are analogous to D-brane or "Wilson bag" excitations in QCD. In Lorentz gauge, the small 
instantons and Wilson line excitations can be described, respectively, in terms of poles and 
cuts of an analytic gauge potential. 



I. INTRODUCTION 



Recent lattice studies of topological charge structure in pure-glue QCD have provided evidence 
that the topological charge density, as constructed from the overlap Dirac operator, organizes 
into low-dimensional long-range structure [jj, consisting of 3-dimensional coherent sheets in 4- 
dimensional spacetime. A further study demonstrated that this structure is inherently global 
and that topological charge fluctuations in QCD should be understood not in terms of individual 
localized lumps but rather as extended brane-like objects P, [J. These lattice results challenge the 
traditional instanton gas or liquid picture Q|, in which the QCD vacuum is described as an ensemble 
of localized (anti-)self-dual lumps of quantized topological charge. On the other hand, the results 
are^uite consistent with the view of QCD theta-dependence provided by string/gauge duality 
, in which topological charge excitations come in the form of membranes or domain walls 
which are the gauge theory dual of wrapped 6-branes in type IIA string theory. Although 
the nature of topological charge excitations in QCD is far from a settled issue, lattice studies 
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may shed light on the dynamics of these coherent fluctuations and their relation to confinement 
and spontaneous chiral symmetry breaking. Analogous long-range sign-coherent 1-dimensional 
regions of topological charge density have also been found recently in a Monte Carlo study of 
2D CP N ~ l sigma models [9|]. Two-dimensional CP N ~ l models and four-dimensional QCD have 
many properties in common. In particular, both are classically scale invariant and have classical 
instanton solutions of arbitrary radius. Lattice evidence for the existence of long-range coherent 
structures of codimension one in both QCD and CP N ^ strongly supports the understanding of 



QCD to^o 
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ogical charge excitations in terms of membranes. In both of these models it has been 
1(1 ] that the observed membranes are associated with boundaries between "k- vacua", 
characterized by an effective local value of 6 which jumps by ±2ir across the boundary. The 
presence of such domain walls in the QCD vacuum was suggested much earlier on the basis of large 
N chiral symmetry arguments [ill ]. The same kind of 3D membrane structure in QCD was also 
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identified by Liischer [12j in the form of "Wilson bags" which are 3-dimensional surface integrals 
over the 3-index Chern-Simons tensor and are analogous to Wilson lines in the CP N ~ l models. 

In Ref. [j| the scaling behavior of the overlap-based topological susceptibility was studied for 
CP 1 , CP 3 , and CP 9 . For CP 3 and CP 9 , xt was found to scale properly with the mass gap. 
However, for CP 1 , xt exhibits anomalous scaling behavior, giving an apparently divergent topo- 
logical susceptibility in the continuum limit. This reproduced a result first obtained by Berg and 
Liischer in a Monte Carlo study of the lattice 0(3) sigma model Q]. The topological suscep- 
tibility xt = {Q 2 ) /V {Q= total topological charge, V= spacetime volume) in two-dimensional 
CP N ~ X sigma models would normally be expected to be a renormalization group invariant param- 
eter, which should scale like the squared mass gap /i 2 (/i = mass of nonsinglet meson) for large 3 
(8 = 1/g 2 , g is the bare coupling constant). However, the Monte Carlo results reported in Ref. 9| 
show that xt does not even approximately scale in the case of CP 1 , although it does scale properly 
for CP 3 and CP 9 . Liischer 14] showed that the anomalous scaling behavior of xt for CP 1 could 



be explained by the contribution of small, point-like instantons with a radius of order a (= lat- 
tice spacing). Although the action of a classical continuum instanton is independent of its radius, 
on the lattice this is not the case, and for CP 1 and CP 2 , the small instantons have such a low 
action that they overwhelm the contribution of more slowly varying fields which might have been 
expected to dominate in the continuum limit. This mechanism is at work for CP 1 and CP 2 , but 
not for CP 3 and higher CP N_1 models, because the action of the small instantons is proportional 
to N and becomes too large to be favored for N > 4. (As we will show, CP 3 is a marginal case 
where instantons are observed but do not dominate as in CP 1 and CP 2 .) The abrupt change in 
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the nature of the dominant topological charge fluctuations as N is increased above 4 provides a 
detailed Monte Carlo example of the phenomenon first discussed by Witten (for both QCD and 



CP N ~ l ) 11], in which instanton effects are exponentially suppressed at large N, and topological 
charge is dominated by quantum effects. This can be described, at least figuratively, by saying 
that instantons "melt" due to quantum fluctuations inherent in the large N c confining vacuum. 
For the case of the CP N_1 models, the instanton "melting point" is N ~ 4. 

The definition of the lattice topological charge density used in the original study by Berg and 
Liischer was an ultralocal discretization of the topological charge operator in the 0(3) sigma model 
formulation. [l^| For our formulation of the lattice CP N ~ l models, the simplest ultralocal definition 
of topological charge is the log of the plaquette constructed from the U(l) gauge links, 

q P (x) = Im{log[^(x)C^(x + /i)tfj(x + /2 + v)Ul(x + z>)]} (1) 

This operator is ultralocal in the gauge fields. It provides an adequate lattice definition of global 
topological charge and gives results for topological susceptibility which are similar but not identical 
to those obtained from the overlap formalism. But as a definition of local topological charge 
density the plaquette-based q p (x) has severe shortcomings. It suffers from a large amount of short- 
range noise associated with local gauge field fluctuations, which obscures the long-range coherent 
topological charge structures found with the overlap-based topological charge density. The overlap 
definition of the lattice topological charge density is Q] 

q(x) = -tx[ l5 D(x,x)\ (2) 



where D is a lattice Dirac operator (in our case the overlap operator |l6j]) which satisfies the 
Ginsparg- Wilson relation 

{D, l5 } = aD l5 D (3) 

This construction of the topological charge density on the lattice provides a more incisive probe into 
the topological structure of gauge field configurations than the plaquette-based operator (^Q). The 
exact lattice chiral symmetry of the overlap Dirac operator leads to an implicit "chiral smoothing" 
of the topological charge which takes place over distances of order a. This reveals the coherent 
structure of the topological fluctuations without any subjective cooling or smoothing procedure. 

In this paper, we investigate the origin of the anomalous scaling behavior of topological suscepti- 
bility in two-dimensional CP 1 and CP 2 sigma models, using the overlap construction of the lattice 
topological charge density. We directly observe small instantons with radii of order a (roughly one 
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to two lattice spacings) in Monte Carlo configurations for both CP 1 and CP 2 . In the continuum 
(large (3) limit these small instantons dominate the topological susceptibility. In contrast to these 
point-like small instanton excitations, the topological charge in CP^ -1 models for N > 4 has been 
found to be dominated by one- d _, M e ^ wan" stetarra a. Afterestablishing 
the prevalence of small instantons in CP 1 and CP 2 , we measure the action of an individual small 
instanton and show that the measured topological susceptibility is in good agreement with a dilute 
instanton gas calculation. Finally we discuss a description of topological fluctuations in Lorentz 
gauge which incorporates both small instantons and domain wall structures into a coherent picture 
in terms of poles and cuts in an analytic gauge potential. 



II. ANOMALOUS SCALING BEHAVIOR OF TOPOLOGICAL SUSCEPTIBILITY 

A. Lattice CP^ 1 Sigma Models 

The Lagrangian for CP N_1 sigma models in the continuum is 

L = d^z-d^i + (z-df.z^iz'jdf.Zj) (4) 

where z l are N complex fields, i = 1, ... ,N, satisfying a constraint z*Zi = 1. This Lagrangian is 
invariant under a local U(l) gauge transformation: i(x), for arbitrary space-time 

dependent a(x). We can introduce a dummy gauge field and rewrite the Lagrangian as 

L = (d ll - iA^z*^ + iAJzt (5) 

where 

A^ = ^(z^z - (0„z)t z ) (6) 

To put CP N ~ l models on the lattice, we introduce 17(1) link fields U(x, x+fi) = e iA ^ x) . CP^ -1 
fields are defined on the sites as z(x). The lattice action consists of gauge-invariant nearest-neighbor 
hopping terms, 

S = (3N ^2 z(x) t C/(x, x + p)z{x + £) + c.c. (7) 

This lattice action is used in our Monte Carlo simulation to generate an ensemble of field config- 
urations. The z fields are updated by a Cabibbo-Marinari heat bath algorithm, while U(l) link 
fields are updated by a multi-hit Metropolis algorithm. 
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B. Overlap Dirac Operator 

The overlap Dirac operator preserves exact chiral symmetry on the lattice via the Ginsparg- 
Wilson relation. For the present work, we employ a standard construction starting from a Wilson- 
Dirac kernel D\y. The overlap Dirac operator can be written as 



L> = ~(l + 75 e(H w (m))) (8) 



where Hw(m) = j^Dw(—m) and 



m i w H w {m) 

e{H w (m)) = (9) 
V H ] w {m)H w (m) 

It can be easily verified that this construction of the overlap Dirac operator satisfies the Ginsparg- 
Wilson relation, Eq. (j3J). 

Using the overlap Dirac operator, we can define the lattice topological charge density by 
Eq. ©, where the trace is summed over spin indices in CP^ 1 models and over spin and color 
indices in QCD. The overlap definition of topological charge is not ultralocal in the gauge links. 
This has the effect of performing some smoothing of the gauge field over a range of order a. For 
smooth U(l) gauge fields in two dimensions, we find that the overlap topological charge at a site 
is reasonably well-approximated by the average of the four plaquette-based charges surrounding 
the site. For a rapidly varying field, and particularly for the small instantons that we observe in 
CP 1 and CP 2 , the overlap-based charge provides a much clearer view of coherent excitations. The 
effect of the smoothing obtained by the overlap construction can also be seen in the structure of 
the 2-point topological charge correlator. Because it is ultralocal, reflection positivity arguments 
imply that the nearest-neighbor correlator for the plaquette-based charge must be negative. On 
the other hand, the overlap-based nearest-neighbor correlator is found to be positive. This allows 
the overlap-based operator to detect sign-coherent structure even if that structure is in the form 
of small 0(a)-radius instantons or membranes of 0(a) thickness. 



C. Topological Susceptibility 

Using the overlap definition of the lattice topological charge density, we have studied the scaling 
behavior of topological susceptibility \t in CP 1 , CP 2 , CP 3 , CP 5 , and CP 9 sigma models. Since 
Xt is generally expected to be a renormalization group invariant parameter in two-dimensional 
CP N ~ l sigma models, it is expected to scale like a (mass) 2 in the continuum limit. There is a 
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natural mass scale in these models which is the mass gap of the z*Zj,(i ^ j) meson correlator 
in the nonsinglet channel. The zero-momentum meson correlator falls off exponentially at large 
Euclidean time, 

J dx2(z*(x)z j (x)z*(0)z i (0)) ~ const, x e"^' 1 (10) 

where u is the mass gap. 

We have measured the mass gap a in lattice units for various values of N and f3 as shown in 
Table 01 We worked on lattice sizes up to 100 x 100. For each N, f3 has been chosen to cover the 



TABLE I: Mass gap of the meson correlator in lattice units 



(3 


CP 1 





CP 2 


& 


CP 3 


(3 


CP 5 





CP 9 


1.0 


.438(5) 


0.9 


.436(5) 


0.8 


.554(2) 


0.7 


.652(6) 


0.7 


.406(2) 


1.1 


.286(5) 


1.0 


.275(2) 


0.9 


.327(3) 


0.8 


.360(3) 


0.8 


.212(2) 


1.2 


.179(3) 


1.1 


.155(1) 


1.0 


.180(1) 


0.9 


.186(3) 


0.9 


.0895(6) 


1.3 


.111(1) 


1.2 


.084(2) 


1.1 


.088(1) 


1.0 


.085(1) 


1.0 


.0579(2) 


1.4 


.070(1) 


1.3 


.043(1) 


1.2 


.0531(3) 


1.1 


.0496(9) 


1.1 


.0475(4) 


1.5 


.036(1) 


1.4 


.0256(8) 


1.3 


.0264(7) 


1.2 


.0308(4) 


1.2 


.0287(4) 


1.6 


.0238(8) 


1.5 


.0207(5) 


1.4 


.0217(5) 


1.3 


.0249(2) 


1.3 


.0255(3) 



region where correlation length £ = l/fi varies from roughly 3 to 50. For each set of N and (3, the 
meson correlator was averaged over 4000 field configurations. The correlator fits were carried out 
by a standard covariant x 2 minimization. 

The topological susceptibility xt can be determined from the large-!/ limit of 

x ( P = (Q 2 )/v (li) 

where Q is the total topological charge on the lattice and V is the space-time volume. An equivalent 
way to measure xt is to integrate the two-point topological charge density correlator, 

x* = J>(*Mo)} (12) 

X 

where q(x) is calculated from the overlap definition (Eq. |2J). If the correlator calculation is carried 
out at finite volume, and the integral is taken over the entire volume of the lattice, this gives results 
identical to those obtained by measuring the global charge Q of each configuration and using Qll|) . 
However, we have found that at large values of (3, where the correlation length is long and finite 
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FIG. 1: The two-point topological charge density correlator for CP 1 with several values of 0. 



volume effects become important, a method based on a truncated integration of the correlator is 
much less affected by finite volume corrections than the calculation based on measurement of the 
global charge Q. The infinite volume topological charge correlator is very short range, consisting 
of a positive core at x < 2 and a short-range negative tail which falls off rapidly, even when the 
nonsinglet meson correlation length is quite long. Fig. ^ shows the point-to-point correlator 
(q(x)q(O)} for CP 1 for several different values of /?. Over this same range, n~ l varies from 2.3 
to 14.4. For larger values of f3, where the correlation length becomes comparable to the box size, 
the global charge Q freezes to zero. The integrated correlator also goes to zero, but, at least for 
correlation lengths up to about half the box size, the shape and magnitude of the short-range 
topological charge correlator remains relatively unchanged from its infinite volume value. The 
fact that the integrated finite volume correlator is smaller comes about through a cancellation 
between the short range part of the integral (which still gives a reasonable approximation to the 
infinite volume susceptibility) and a negative background contribution which does not fall off as 
rapidly with distance and therefore gets most of it's contribution from larger \x\. By comparing the 
correlator on several box sizes, we find that this longer range negative background is a finite volume 
effect. As a result, on smaller boxes the integrated short-range correlator gives a better estimate 
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of the infinite volume xt than one could obtain from the fluctuations of the global charge. After 
experimenting with different prescriptions for cutting off the short-range correlator integration, we 
have chosen to calculate the correlator as a function of \x\ and integrate out to a cutoff value of 
\x\ < x c for which the measured correlator is within one standard deviation of zero. This was 
always in the range 3 < x c < 5. Results were essentially unchanged if we used a constant cutoff of 
x c = 4. To test the method, we calculated xt as a function of the cutoff radius for CP S at /3 = 1.0 
on various lattice sizes, as shown in Fig. [3 Xt measured on the 40 x 40 lattice matches nicely with 
that measured on the 50 x 50 lattice at each cutoff radius. This shows that the 50 x 50 lattice is 
large enough to give a correct xt that is close to the infinite volume susceptibility. By using (|TT|) 
on the 50 x 50 lattice configurations, we obtained xt = 0.0032(3). From the graph, we observe that 
Xt measured on the 40 x 40 or 50 x 50 lattice quickly converges to the infinite volume result when 
x c ~ 4. However, on the 20 x 20 lattice, xt converges to a value lower than the infinite volume 
result. On the other hand, if we apply the cutoff at x c ~ 4, we can still extract approximately the 
right value of xt on the 20 x 20 lattice. Most of the finite volume effect that causes the globally 
determined xt to be too small on the 20 x 20 lattice comes from distances > 4. On the larger 40 x 40 
and 50 x 50 lattices, the value obtained for xt is essentially independent of the cutoff for x c > 4, 
indicating that the longer range negative contribution from |x| > 4 on the smaller box is a finite 
volume effect. This shows that the truncated integration method is less affected by finite volume 
than a calculation based on measurement of the global charge Q. Using the truncated integration 
method, we have measured xt in lattice units for various values of N and (3, as shown in Tabic ITT1 



TABLE II: Topological susceptibility \t in lattice units 





CP 1 


& 


CP 2 


a 


CP 3 


P 


CP 5 


a 


CP 9 


1.0 


.0146(2) 


0.9 


.0137(2) 


0.8 


.0147(2) 


0.7 


.0154(3) 


0.7 


.0051(1) 


1.1 


.0106(2) 


1.0 


.0081(2) 


0.9 


.0079(2) 


0.8 


.0069(1) 


0.8 


.00105(5) 


1.2 


.0070(2) 


1.1 


.0040(1) 


1.0 


.0034(1) 


0.9 


.0021(1) 


0.9 


.00021(2) 


1.3 


.0040(1) 


1.2 


.00157(7) 


1.05 


.00194(8) 


1.0 


.00027(4) 


1.0 


.00009(2) 


1.4 


.0022(1) 


1.3 


.00056(5) 


1.1 


.00100(6) 


1.1 


N/A 


1.1 


N/A 


1.5 


.00104(5) 


1.4 


.00035(3) 


1.15 


.00057(5) 


1.2 


N/A 


1.2 


N/A 


1.6 


.00053(5) 


1.5 


N/A 


1.2 


.00026(4) 


1.3 


N/A 


1.3 


N/A 



To show the scaling behavior of xt with respect to (mass) 2 , we plot Xt/^ 2 vs. // in Fig. A 
constant value on the graph indicates proper scaling. From the graph we observe that xt appears to 
scale properly for CP 3 , CP 5 , and CP 9 in the continuum limit, but it does not even approximately 
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FIG. 2: xt vs. cutoff radius 




FIG. 3: Scaling behavior of Xt/y 2 in CP 1 , CP 2 , CP 3 , CP 5 and CP 9 
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scale for CP 1 and CP 2 . Our results confirm that the overlap-based xt exhibits the same anomalous 
scaling behavior observed in the 0(3) sigma model formulation of CP 1 |l^ |. From Fig. ISlwe see 
that this anomalous scaling behavior occurs in CP 1 and CP 2 , but not in CP 3 , CP 5 , or CP 9 . This 

n 

is exactly the behavior predicted by Liischer who showed that, for N < 3, the topological 
susceptibility should be dominated by small instantons with a radius of order a. In the next section, 
we will discuss the theoretical explanation of the anomalous scaling behavior of xt, present the 
direct observation of small instantons in Monte Carlo configurations, and discuss some properties 
of these objects. 



III. SMALL INSTANTONS 

The divergent behavior of xt m CP 1 and CP 2 can be attributed to instantons with a radius 
of approximately one to two lattice spacings. These small instantons are short-range gauge field 
fluctuations with quantized topological charge. In the space of lattice gauge field configurations 
they reside near the boundary between the Q = 1 and Q = sectors, in the sense that a small 
adjustment of the gauge links will cause the instanton to "fall through the lattice" and turn into 
a Q = configuration. The dominance of these small instanton configurations violates the usual 
assumptions about the continuum limit. In a naive view of the continuum limit, one would assume 
that configurations are dominated by fields for which the values of gauge-invariant local densities 
(e.g. action density or topological charge density) are smooth over distances of order a. Although 
the probability for the small instantons to occur is exponentially small, the probability for more 
smoothly varying fields to have nonzero total topological charge is also exponentially small. Which 
type of fluctuation dominates the topological susceptibility depends on the precise value of the 
instanton action. 

If we only take into account the contribution from the slowly varying fields, the topological 
susceptibility should scale like a (mass) 2 in the continuum limit, i.e. 

X r oc (3*/ N e-W, (P^oc) (13) 

However, if we take into account the contribution from small instantons, a dilute gas calculation 
yields 

X^-oc/rV*, (0-oo) (14) 
where (3e is the minimal action of gauge field configurations in the Q = 1 sector. 
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In Ref. ^| the minimum-action field configuration with Q = 1 was constructed and its action 
was calculated. This gave the value of e in the 0(3) model as e = 6.69 . . . For other CP^" 1 models, 
it was argued that the minimal action of a small instanton should be proportional to N, according 
to the inherent factor of N in the CP N ~ l action (see Eq. [7J). Since e < 4-7T for CP 1 and CP 2 , the 
small instantons make the dominant contribution to the topological susceptibility. This explains 
why xt does not scale as /x 2 in the continuum limit in CP 1 and CP 2 . However, for CP^ -1 models 
with N > 3, e = y • 6.69 > An, so xt is dominated by less singular topological fluctuations. On the 
other hand, even for N > 3, the dominant excitations are not smooth distributions of topological 
charge, but rather, singular domain-wall type excitations which give rise to a dominant contact 
term in the 2-point correlator 

Even a qualitative survey of overlap topological charge distributions in the Monte Carlo config- 
urations is sufficient to identify small instantons as the origin of the divergence of xt m CP 1 and 
CP 2 . To emphasize this point, we will first present a few typical configurations. Small instantons 
are particularly prominent and easy to identify at large values of (3, where the meson correlation 
length is large. Interestingly, the finite volume effects on the fluctuation of the global topo- 
logical charge Q are completely different for CP 1 and CP 2 than they are for the N > 4 models. 
For CP 3 and higher, the fluctuation of Q to nonzero integer values "freezes out" at large values 
of (3 for a lattice with a fixed number of sites. When the correlation length becomes comparable 
to the box size, Q stops fluctuating altogether. In marked contrast, for example, CP 1 on a 40x40 
lattice continues to fluctuate to nonzero Q well above (3 = 1.6, where the correlation length 
is greater than 40. This is apparently due to the fact that the small instantons remain of order 
lattice spacing in size, even when the correlation length becomes large. They are therefore not very 
susceptible to being squeezed out by finite volume effects at the scale of the correlation length. 



In Fig. 0] we show several typical configurations. Fig. 4(a) is a CP 1 configuration with (3 = 2.5 
and Q = 1.0. It clearly shows a positively charged peak which is distinctly different from the 
background quantum fluctuations. We will refer to these objects, which are easily identified in 
the large (3 regime, as "peak structures." As we discuss below, the amount of charge contained 
in a typical peak structure is approximately one unit, so they clearly exhibit local quantization 
of topological charge. Furthermore, for large values of j3, fluctuations of the global charge Q by 
a unit during the Monte Carlo run are always accompanied by the appearance or disappearance 
of one of these peak structures. For comparison, we also show a CP 1 configuration with (3 = 2.5 
and Q = in Fig. |4(b)[ where no peak structure can be seen. Fig. 4(c) shows two negative peak 



structures in a CP 1 configuration with f3 = 1.6 and Q = —2.0. Because of the smaller value of (3, 
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(e) CP 3 , 13 = 1.6, Q = (f) CP 9 , f3 = 0.9, Q = -1.0 

FIG. 4: (a) A peak structure is clearly seen in a CP 1 configuration with (3 = 2.5 and Q = 1.0. (b) No peak 
structures in CP 1 with (3 = 2.5 and Q = 0. (c) Two negative peak structures in CP 1 with (3 — 1.6 and 
Q = -2.0. (d) A peak structure in CP 2 with (3 = 1.8 and Q = 1.0. (e) CP 3 with /3 = 1.6 and Q = shows 
no peak structure, (f) CP 9 with /3 = 0.9 and Q = —1.0 shows no peak structure. 
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the background quantum fluctuations are seen to be larger than in the previous plots. The CP 2 
configuration in Fig. 4(d) also shows a peak structure with (5 = 1.8 and Q = 1.0. However, at 
comparable correlation lengths, the CP 3 and CP 9 configurations contain no prominent peaks and 



always have zero total topological charge on a 40 x 40 lattice. Fig. 4(e) shows a CP 3 configuration 
with j3 = 1.6 and Q = 0, where no peak structure can be seen. Finally, Fig. |4(f )| shows a CP 9 
configuration with f3 = 0.9 and Q = —1.0. Here, in spite of the fact that Q is nonzero, there is 
no sign of an integer charged peak structure. In fact, for CP 5 and CP 9 , there were no integer 
charged peak structures in any of the configurations studied. The case of CP 3 is interesting be- 
cause, although there is some indication of small instantons (see below), they are not as dominant 
as in CP 1 and CP 2 . Note that the one-instanton action of e = 14.1 for CP 3 is slightly above the 
instanton "melting point" of 4n = 12.57. . . predicted in the dilute gas approximation. In the case 
of CP 3 , the small instantons do not lead to any observable anomalous scaling of \t because their 
action is very close to 47r. 

To add further support for the assertion that the peak structures in CP 1 and CP 2 are indeed 
small instantons, we measured the topological charge inside the highest peak structure in each 
field configuration. As we will show in the following discussion, a good effective definition of 
the topological charge within a peak is given by the sum of the charges within a radius < 2 of 
the maximum-charge site. To illustrate the quantization of the charge in the peak structures, 
we analyze, for several values of (3 and N, a large sample of configurations which have a global 
topological charge of Q = ±1. For each configuration we calculate the amount of charge within the 
highest peak structure. The histograms show the number of configurations which have a charge 
within a given range for the highest structure. As clearly shown, the CP 5 and CP 9 models exhibit 
very different histograms from CP 1 , CP 2 and CP 3 . The CP 5 and CP 9 configurations typically 
have much less than a unit of topological charge in the highest peak structure, while in CP 1 , CP 2 
and CP 3 , the majority of Q = ±1 configurations have a structure which contains approximately 
a unit of topological charge. This clearly confirms that the peak structures we observed in Fig. 0] 
are in fact small instantons with a radius of < 2 lattice spacings. By contrast, the topological 
structures in higher CP^ -1 models exhibit no tendency toward local quantization, consistent with 
expectations from the domain-wall picture of topological charge excitations in these models. 

To estimate the approximate radial size of the peak structures, we considered an ensemble of 
CP 1 configurations on a 40x40 lattice at f3 = 1.8, where we have 72 Q = ±1 configurations out 
of 400. In each configuration, we started with the highest peak, then measured the topological 



charge density q(r) and the integrated topological charge Q(r) within radius r. Fig. 6(a) shows 
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FIG. 5: Histogram of the number of configurations with respect to the topological charge inside the highest 
peak structure 



the r dependence of q{r) and Q{r) with error bars. We see from the graph that q(r) drops to 
zero at about r ~ 2, while at the same time Q{r) goes up to ~ 1.0. For this reason, we define 
the peak structure as occupying the lattice sites within a radius of two lattice spacings from the 
peak. We also observe from the graph that the q(r) fluctuates around zero when r > 2.0 and the 
fluctuation is comparable to error bars. This confirms that the radial size of the peak structures 
we observed in Fig. 0] is indeed roughly 1 to 2 lattice spacings. All of these properties fit precisely 
the description of small instantons, which arise on the lattice as minimum action configurations in 
the Q = ±1 sector. 

In order to compare topological susceptibility measurements with the predictions of a dilute gas 
of small instantons, we want to determine the action of a single small instanton by measuring the 
minimal action in the Q = 1 sector for various CP N ~ 1 models. The value of the minimal action 
for a small instanton determines whether they will dominate the topological fluctuations in the 
continuum limit. For CP 1 and CP 2 the value of the minimal action also determines the form of 
the divergence of Xt m the continuum limit. We constructed a small instanton configuration on 
the lattice, which consists of three vertical links of e lT1 , e ln , and e~ lTI respectively on three adjacent 
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FIG. 6: (a) The r dependence of the topological charge density q(r) and the total topological charge Q(r) 
within radius r; (b) The gauge field action steps down when the CP 1 configuration tunnels from Q = 1 into 
Q = sector. The upper curve is the action and the lower curve is the charge Q. 



sites on a horizontal row, while all other links are set to one. Here the value of r] is taken at the 
boundary between the Q = 1 and Q = sectors, as determined by the overlap Q. This gives, 
rj = .0873(1). We then annealed this configuration (by running the Monte Carlo on the z fields 
at large j3 with fixed gauge links) and measured its action. The minimal actions we obtained are: 
e = 6.76 for CP 1 , e = 10.41 for CP 2 and e = 14.08 for CP 3 . For CP 1 this is close to Liischer's 



result e = 6.69 [Tj]. Since our lattice action and definition of topological charge differ from the 0(3) 
model formulation, exact agreement between our measured e and that of Ref. 14] is not expected. 
We also looked at how the lattice action changes when we let a CP 1 configuration tunnel from the 
Q = 1 sector into the Q = sector, starting with a Monte Carlo generated Q = 1 configuration 
and cooling it by running at a very large j3. An action step is expected which should be close to 
the minimal action in the Q = 1 sector. Fig. |6(bJ| clearly shows the action step when the tunneling 
occurs. The action step is measured to be about 7, which is very close to the one-instanton action. 



Some insight into the dynamics of small instantons in CP 1 and CP 2 can be obtained by compar- 
ing the measured topological susceptibility with the prediction of a dilute instanton gas calculation. 
In the dilute gas approximation, the topological susceptibility is determined by the one-instanton 
action e, and behaves like 

xf 9 - oc /TV* 09 - oo) (15) 
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FIG. 7: Plot of the coefficient 8 t {(3) denned in Eqn.EO 
To test how well the dilute gas approximation works, we plot the quantity 



(16) 



against (5 in Fig. [7[ up to an irrelevant overall factor. Fig. 7(a) and |7(b)| show that 5 t becomes 
roughly constant at large (3 for CP 1 and CP 2 , confirming the dilute gas calculation. However, 5t 
is a rapidly increasing function of (5 for CP 5 and CP 9 , showing that the expected contribution 
from instantons (|15|) is much smaller than the measured topological susceptibility. This confirms 
the conclusion that small instantons are irrelevant to the determination of xt for N > 4. 

In the next section we discuss the mechanism of topological charge fluctuations in different 
(jpN-i moc i e i S- \y e w [\\ ar g U e that small instantons and membrane structures are the two types of 
lower dimensional topological structures (zero- and one-dimensional respectively) that can appear 



17 



in CP N 1 models. 

IV. MECHANISM OF TOPOLOGICAL CHARGE FLUCTUATIONS 

As we have shown in the previous section, small instantons dominate topological charge fluc- 
tuations in CP 1 and CP 2 and produce anomalous scaling behavior of xti while CP 3 and higher 
QpN-l moc i e i s exhibit proper scaling behavior xt ~ ^ 2 ■ Thus, in the continuum limit, CP 1 and 
CP 2 have infinite topological susceptibility. The original observation of this divergent susceptibil- 
ity was obtained with an ultralocal definition of topological charge in the 0(3) sigma model 
formulation. The fact that the same divergent behavior of Xt/ 1^ 2 occurs when one uses the overlap 
construction suggests that the divergence cannot be eliminated by a better lattice definition of 
topological charge, and that it has some physical significance in the continuum limit. To gain a 
better understanding of this, it is useful to consider the topological charge correlator in Euclidean 
2-space, 

G(x) = (q(x)q(0)) (17) 

The lattice studies discussed here and in Ref. [jj provide clear evidence that, in the continuum 
limit, this correlator is very short range (compared to the correlation length /i _1 ), so one might 
expect an operator product expansion to provide insight. Here we consider an operator product 
expansion in the effective U(l) gauge theory obtained by integrating out the z-fields. Taking the 
large N solution as a guide, we expect that the charged z-particles acquire a "constituent quark 
mass" M. At the same time, z-loop effects produce a dynamically generated kinetic F 2 term in 
the effective gauge action. This effect leads to a confining Coulomb potential between constituent 
z-particles. However, in the large iV limit, the confining potential is weak, corresponding to a 
string tension of order M 2 /N. Thus, one can distinguish between a relatively long distance scale 
~ associated with confinement, and a shorter distance scale of order M _1 , which probes 

the substructure of the constituent z particles. By invoking the operator product expansion in the 
effective gauge theory obtained by integrating out the z-fields, we are assuming that the distance 
scales relevant to topological structure are characterized by the confinement scale. The large iV 
solution, discussed below, provides a good illustration of the order of limits involved here. In the 
large ./V approximation, after integrating out the z fields, the action L(x) can be replaced by 

L(x) - \F% (18) 
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where Fn U is a rescaled field strength obtained by the replacement 



UttM 2 a . . 

^^V~ A " (19) 

The two lowest dimension gauge invariant, Lorentz invariant local operators that can appear in 
the OPE for (|17f) are the identity and the squared field strength F 2 = F 2 , 

q(x)q(0) ~ d(x) + C 2 (x)F 2 (0) + ... (20) 

where . . . represents operators of dimension greater than 2. In the following discussion, it will 
also be useful to consider the corresponding OPE coefficients for the product of two Chern-Simons 
currents. Defining 

tf S (?) = ^uA^(x) (21) 
we write the OPE (ignoring possible gauge-dependent terms which do not contribute to (|2UlO . 

S (x)jC s (0) ~ Ci^(x) + C 2flu (x)F 2 (0) + ... (22) 
The Chern-Simons correlator 

g™{x) = (J c ; s {x)jc s m (23) 

is related to the topological charge correlator (|TT)) by 

G(x) = -d^G^ix) (24) 

In the free theory defined by (|18j) . the behavior of the OPE coefficients is given by naive dimensional 
counting. For the topological charge OPE, this would give. 

d(x) ~ 1/x 4 (25) 
C 2 {x) ~ 1/x 2 (26) 

However, we will see in the following discussion that, although such power-law behavior is true 
in the OPE for the Chern-Simons correlator, the corresponding OPE for the topological charge 
correlator "collapses" to pure contact terms of the form 

Ci{x) ~ V 2 5 2 (x) (27) 
C 2 {x) ~ 5 2 {x) (28) 
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As we argue below, this collapsing of the operator product expansion for G(x) is a consequence of 
the lower dimensional character of the small instanton and domain wall excitations of topological 
charge. The Fourier transformed OPE coefficients 

Ci(q) = J tPxdWe** (29) 

have large q 2 behavior 

Ci{q 2 )~q % (30) 

C 2 (q 2 ) ~ 1 (31) 

The singular terms in the operator product expansion correspond to polynomial subtraction terms 
in the dispersion integral representation of the two-point correlator. We consider the correlator in 
momentum space, 

G(q 2 ) = [ d 2 xe iqx G{x) (32) 



Under general assumptions, the correlator can be written in a dispersion integral representation 

3, 
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i2 P(A 2 ) 



G(q 2 ) = Cl q 2 + c 2 - / d\ 2 -^—^ (33) 
Jul T + * 2 



(Note: This representation holds in two space-time dimensions. In four dimensions, the polynomial 
part can include a q 4 term.) The first two terms represent contact terms associated with the short 
distance singularities of the OPE 1)20(1 . The integral term represents the contribution of real 
intermediate states of invariant mass A. Here M$ is the mass gap in the flavor singlet channel. 
The standard way of analyzing correlators on the lattice is to pick a time direction and Fourier 
transform over the spatial direction with spatial momentum component q\ = 0, i.e. 

G(t) = J d Xl G(x = t, Xl ) (34) 

The dispersion representation then gives 

poo j\2 

G(t) = - Cl 5" (t) + c 2 5(t) - / p (A 2 ) e - A l*l (35) 

JMl 2A 

The topological susceptibility is 

poo J\2 

X ' = C2 -Jm S ^ P(X2) (36) 
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Note that the c\ contact term arises from the identity operator term in the OPE, while the C2 
term arises from the F 2 insertion term. The c\ term does not contribute to Xt an d the C2 term 
must be positive and larger in magnitude than the negative contribution from the integral over 
real intermediate states in order for \t to be positive. 

A detailed numerical analysis of the spectral and scaling properties of the q(x) correlator will 
be reported elsewhere. Here we will quote and utilize some of the main properties of the correlator 
which emerge from this numerical study. The results are found to be consistent with the assumption 
that the correlator is reasonably well-described by contact terms alone, without the need for an 
intermediate state term. We should emphasize that a description of the correlator as purely contact 
terms does not mean that real propagating intermediate states do not contribute, but only that 
those states are heavy compared to the q 2 values we are probing on the lattice when we calculate 
the two-point correlator. The effect of heavy intermediate states can be absorbed into the contact 
terms via 

\2 



C-2 



c i + / -rrP (A 



C-2 



The large N solution 
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19] 



Ml 



"A 2 " 



(37) 
(38) 



provides an explicit example which illustrates how the pure contact 
term approximation to G arises. We consider the gauge field vacuum polarization tensor (the inverse 
of the (ApAJ) correlator) in momentum space, 



A 



[IV 



Since 2i:q(x) is just the curl of A 



q(x) 



1 

2^ 



£uvdn,A u 



the gauge correlator is related to the topological charge correlator by 

,2 



(40) 



(41) 



or 



4vr 2 G(q 2 ) = q 2 A(q 2 ) (42) 
In the large N limit, this is determined by a one loop calculation depicted in Fig. |H1 which gives 



A(q 2 



2vr 

iV 



e(9 2 )In 



C(Q 2 ) - 1 



(43) 
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FIG. 8: One-loop graphs contributing to the gauge field correlator in the large N approximation. 



where £,(q 2 ) = y 1 + ^p-- In this example, the pure contact term approximation to the topological 
charge correlator is obtained by assuming that the mass gap Mq = 2M is large compared to the 
momenta relevant to topological charge structure, i.e. AM 2 » q 2 . Expanding in powers of q 2 /M 2 
and neglecting terms that vanish in the limit q 2 /M 2 — ► 0, we get 

Note the presence of a q 2 = pole in the gauge field correlator. Since the Chern-Simons current is 
just the dual of the gauge field, the same massless pole also appears in the Chern-Simons correlator, 
and it's residue determines the topological susceptibility in the large iV limit to be 

*-!£ + o(l) (45, 



irN \N 2 , 

When formulated in terms of the Chern-Simons correlator, 2D CP N_1 and 4D QCD have a similar 
topological charge structure at low momentum. In both cases, the topological susceptibility is 
just the residue of the massless pole in the Chern-Simons current correlator. The massless Chern- 



121. 



Simons pole embodies the "secret long range order" associated with topology in gauge theory 
The physical significance of this massless pole is subtle. Since the CS correlator is not gauge 
invariant, a q 2 = pole does not imply the existence of a massless particle. On the other hand, 
since the residue of the pole is the gauge invariant topological susceptibility, the pole term in the 
CS current correlator is in fact gauge invariant and cannot be transformed away. It thus represents 
a real physical long range correlation which is built into the gauge field by virtue of its topological 
fluctuations. Nevertheless, in the q 2 « AM 2 approximation, the topological charge correlator 

G(q 2 ) = ^ 2 A(g 2 ) (46) 

is a polynomial in q 2 , so that in coordinate space it collapses to local contact terms. Numerically 
on the lattice, this leads to a G{x) that is extremely short-range. 
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In the remainder of this section, we will discuss the contribution of small instantons to the 
correlator in the approximation of keeping only the contact terms in G(t). In this approximation, 
the time-dependence of the correlator is simply 



G(t) = - Cl 6 (t) + c 2 6(t) 



(47) 



We note that the numerical calculation of the topological susceptibility from the integral of the 
lattice correlator 9] involves a large cancellation between the positive part of the correlator at t = 
and 1, and the negative short-range tail from t = 2toi~4or5. A similar cancellation has been 
observed in 4D QCD 3]. For the CP N ~ l models, we find that, as beta is increased, the range of 
the negative tail remains approximately fixed in lattice units as discussed in Section II. This leads 
us to conclude that the negative tail of the correlator may be treated as part of the contact term, 
arising from the 5 term in (|47|). 

Normal scaling behavior for the coefficients c\ and C2 would be 



c\ ~ const. 



C2 ~ /i 



(48) 
(49) 
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In order to use Eq. (|47|) as a fitting formula to extract the coefficients of the contact terms from 
the Monte Carlo data, we employ the following parametrization of a smeared delta function: 

'<*> - m 

The second derivative of this expression also provides a fitting function for the S"(t) term in the 
correlator. We then fit the correlator to a sum of a 5(t) and a S"(t) term. The parameters which 
are allowed to vary in these fits are c\, the coefficient of the 5" term, and d, the smearing distance. 
The parameter C2 is just the susceptibility Xti which is computed separately and held fixed in the 
fits. With these lattice approximations to 5(t) and S"(t), the expression (|47|) provides a good fit 
to all of the correlators studied. Fig. shows a typical fit. This plot is for CP 1 at (3 = 1.6. The 
scaling behavior of the smearing parameter d provides a significant test of the pure contact term 
assumption, Eq. (|47[). In Fig. ^]we plot the fit values for d as a function of n for various CP N ~ 1 
models. In all cases, d approaches a finite value in lattice units of d ~ 1.5 lattice spacings. This 
is roughly the same as the measured thickness of domain walls reported in Thus, the physical 
range of the correlator goes to zero in the continuum limit. This provides strong numerical support 
for the assumption that the correlator is well-approximated by contact terms alone. The scaling 
properties of C2 = Xt for the various CP N_1 models have already been discussed. In particular, 
for CP 1 and CP 2 , C2 exhibits anomalous scaling properties due to small instantons. In Fig. Illl we 
have plotted the value of the c\ coefficient as a function of \i for several CP N_1 models. From these 
results, we conclude that c\ scales properly (i.e. approaches constant) as a function of \i for all the 
CP models. Most notably, unlike C2, the coefficient c\ does not appear to receive divergent 
contributions from small instantons, even for CP 1 and CP 2 . The result that small instantons 
contribute only to the C2 coefficient and not to the c\ coefficient is also expected on theoretical 
grounds, as we discuss in the next section. The fact that c\ scales canonically (i.e. approaches 
a finite constant) in the continuum limit provides further support for the consistency of the pure 
contact term approximation for the topological charge correlator. 



V. HOLOGRAPHIC AND CONFORMAL PROPERTIES OF TOPOLOGICAL CHARGE 

IN CP*- 1 

All of these results have a natural interpretation in terms of a two-component picture for topo- 
logical charge fluctuations in CP^ 1 models, which we propose and discuss in this Section. To 
motivate this picture, we first observe that the required negativity of the (q(x)q(O)) correlator for 
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nonzero separation has a fundamental dynamical effect which strongly constrains the types of topo- 
logical charge fluctuations which can contribute significantly to the vacuum path integral. The way 
that this constraint is realized in Monte Carlo configurations is by the predominance of topological 
charge distributions which have a "subdimensional" character, i.e. ones in which coherent regions 
of topological charge are effectively either zero-dimensional (small instantons) or one-dimensional 
(Wilson lines). A more unified perspective on both of these types of topological charge excitations 
can be obtained by considering Witten's holographic formulation of theta dependence in gauge 
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theory In this formulation, the theta term in 4D Yang- Mills arises from compactification of a 
5D Chern-Simons term. Similarly, we can interpret the theta term in CP N_1 as a compactified 
3D Chern-Simons term, 

Cos = ie abc A a d b A c (51) 

Here, a,b,c, . . . run from 1 to 3. Let us denote the original spacetime dimensions by 1 and 2, 
and the compactified dimension by 3. Then in the limit of small radius of compactification, the 
Chern-Simons term reduces to a theta term, 

Cos -> = i9q{x) (52) 

where fi, v = 1, 2, and 

' A 3 dx 3 (53) 



From this 3-dimensional framework, a small instanton in CP 1 or CP 2 can be interpreted as a 
charged particle coupled to the gauge field ^3 which has a world line wrapped around the compact 
direction in a closed loop, and is pointlike in the 1-2 plane. On the other hand, we may integrate 
by parts and write, 

Cos = —e^d^Av = J V A V (54) 

where 

Jv = -r-e^dnfi (55) 

In this way of writing the CS term, the current couples to the gauge field in the 1-2 plane. In 
the limit of small compactification radius, the quantity 9 defined by ()53j) reduces to the constant 
theta parameter of the 2D theory, but only mod 2irk, where k is an integer which labels a local 
k-vacuum state. These k-vacua are separated by domain walls, and the current is an "edge 
current" which is nonvanishing along these domain walls. 

Thus, the 3D Chern-Simons/holographic view of theta dependence in CP N ~ l models provides 
a more unified picture of the two dominant types of topological charge excitations in CP N ~ l 
models. In 3-dimensional spacetime, both small instantons and the domain walls associated with 
Wilson line excitations are 0-branes, i.e. particle-like excitations with 1-dimensional world lines. 
They differ only in the orientation of their world line with respect to the compactified dimension, 
with small instantons being wrapped around the compact direction localized in the 1-2 plane, and 
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domain walls being world lines stretched out in some direction in the 1-2 plane. Some further 
insight into this structure, which also exposes an apparent string/gauge connection, can be gained 
by considering the 2D CP N_1 models in Lorentz gauge, d^A^ = 0. In this gauge, the 2D field 
may be written in terms of a single scalar field: 

A„ = (56) 

The Chern-Simons current is the gradient of and the topological charge is its Laplacian, 

V 2 $(x) = 2irq{x) (57) 

The Monte Carlo results discussed here and in Ref. ($]) suggest the following idealization: Let us 
suppose that regions of nonzero topological charge in the vacuum are confined to zero- and one- 
dimensional submanifolds of two-dimensional spacetime. In this idealized view, spacetime consists 
mostly of voids in which q{x) = 0, with topological charge being confined to either isolated points 
(small instantons) or to boundaries between voids (domain walls). In the voids q{x) = 0, and (|57|) 
implies that &(xi,X2), within a localized region of spacetime and for a given gauge configuration, is 
a harmonic function, which can be written in terms of holomorphic and anti-holomorphic conformal 
fields, 

$( Xl ,x 2 ) = <t>{z) + <t>{z) (58) 

where z = x\+ix2- In this framework, conformal symmetry is broken by the presence of instantons 
and domain walls in the vacuum, represented respectively by poles and cuts in the analytic structure 
of 4>(z) configurations contributing to the path integral. The local conformal properties of the gauge 
field implied by this picture provide some insight into the collapsing of the OPE for the topological 
charge correlator discussed in Section IV. To illustrate this point, consider the leading OPE term 
for the CS correlator (|22[). which is conformally invariant. This vacuum insertion term should be 
equal to the exact correlator in the absence of conformal symmetry breaking. The second term 
in (|22|) is proportional to the (dimensionful) topological susceptibility, and thus arises from the 
breaking of scale invariance due to the presence of poles and cuts in the analytic structure of cp(z). 
In terms of the Lorentz gauge scalar field $, the Chern-Simons current is J^ s = d^&, and the OPE 
coefficient C\^ u {x) is completely determined up to an overall constant by conformal invariance of 
the holomorphic and anti-holomorphic terms in the correlator, 
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and a similar expression for the anti-holomorphic part. The cross-correlators between a holomor- 
phic and an anti-holomorphic field vanish. Using (I56j) and (|58|h (|59j) and the analogous correlator 
for dz<j)(z) completely determine the form of the leading OPE coefficient for the Chern-Simons 
correlator up to an overall constant, 

CV(s) = const. x ^-yV (6Q) 

It can now be seen that this OPE coefficient collapses upon differentiation, i.e. 

d t j,Ci llu (x) = const, x d v 5 2 {x) (61) 

This vacuum insertion term produces the 5" (i) term in the G(t) correlator but does not contribute 
to xt- 

In the topological charge correlator, the F 2 insertion term (second term in the OPE (|2U|) or 
l|22jl) breaks conformal invariance and provides nonvanishing xt- Again, dimensional considerations 
and Lorentz invariance allow us to write 

CV = const, x Q<5 M „ lnx 2 + (62) 

(Here we have dropped an infrared divergent constant oc S^lnX 2 , where A is an infrared cutoff 
mass.) The relative coefficient of the two terms in (j62j) is determined by the Lorentz gauge condition 
on the dual correlator. The long range lnx 2 behavior of the Chern-Simons correlator coming 
from (|62j) corresponds to the massless 1/q 2 pole in the momentum space correlator (|43|) or (|44|). 
Thus, the large x behavior of (|62j) exhibits the hidden long range order of the theory. Again, the 
corresponding term in the topological charge correlator collapses to a contact term 

C 2 (x) = -d^d u C 2 ^(x) = const, x 5 2 {x) (63) 

It is this term in the correlator which gives nonzero topological susceptibility. 

In Section IV it was shown numerically that small instantons only contribute to the C 2 coef- 
ficient and not to the C\ coefficient in the OPE for the q{x) correlator. We can now present a 
theoretical argument for this by showing that the contribution of small instantons to the Chern- 
Simons correlator is of precisely the form (|62j) . The distribution of Chern-Simons current around 
a small instanton at spacetime point y is given by 

J cs (x) - 1 ( x ~ y ^ (qa) 
J » [ j "2vr(x-y)2 [M) 
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Thus, the contribution of an instanton to the correlator ( J^' s (x/2) J„ s (— x/2)) is proportional to 

d 2 y {y -\x)^{y + \x), 



- h + h (65) 

where 



(27T) 2 (y- ^(y+l^ 



/ 



(2vr) 2 (y-x/2) 2 (y + x/2) 2 
and 



(66) 



1 f d 2 y 1 

h -~~^ Xv JWf {y-x/2)\y + x/2f (67) 

These integrals are of the same form as the one-loop Feynman integrals for the large- iV correlator 
(although here the integral is over instanton moduli space rather than momentum space). Again, 
as in (|62j) . we drop infrared divergent constants which do not contribute to the topological charge 
correlator. This gives 



7 "i(^-^) (lnx2 - 2) m ' 



and 



l2= lx^ lnx 2 (69) 
8-7T X Z 



Combining these (and dropping a constant term oc 5^) we get 



h+h = -iz(^- 2 + J S L ) (70) 



Thus a small instanton gives a contribution to the correlator of precisely the same form as the 
coefficient, (|62j) . We conclude that, within the dilute gas approximation, small instantons 
contribute only to the C2 coefficient, in agreement with the numerical results of Section IV. 



VI. CONCLUSIONS 



Monte Carlo calculations using the overlap Dirac operator to study topological charge structure 
in CP N ~ l models have begun to clarify some longstanding issues in these models. The CP N ^ 
models offer an ideal laboratory for studying topological charge in asymptotically free gauge theory, 
and precise and detailed parallels can be drawn between these models and QCD. Witten's original 

n 

arguments that instantons in QCD will "melt" due to quantum fluctuations drew heavily on 
the CP N ^ analogy, in particular the absence of instanton effects in the large iV solution and 
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nonzero topological susceptibility at order 1/N. Liischer's arguments that topological charge in 
QCD could be associated with "Wilson bag" surfaces were also supported by the CP N ~ l analogy, 
where the corresponding excitation is just the familiar Wilson line. Perhaps the most important 
simplification in the CP N ^ models, relative to QCD, is the fact that in one spatial dimension, 
a domain wall is a point-like excitation. This allows the success of standard field theory methods 
such as the large N expansion. The world lines of constituent z-particles which emerge in the large 
N limit can be identified with domain walls between k- vacua, where the value of 0/2tt, interpreted 
as an electric field, changes by one unit of flux. 

For the CP N ~ l models, a value of N « 4 marks the instanton melting point, i.e. the transition 
from dominance of q{x) by O-dimensional instantons to dominance by one-dimensional line-like 
excitations. From a holographic 3D Chern-Simons viewpoint, the melting of an instanton is actually 
the unwinding of a Wilson line excitation which was wrapped around the compact dimension. 
Thus the line-like excitations which are found to dominate the larger N models might be viewed 
as not only a replacement for the instantons, but in some sense the remnants of melted (unwound) 
instantons. 

The full dynamical implications of small instantons in the CP 1 and CP 2 models remain to be 
explored. One immediate question is whether these models can have infinite topological suscep- 
tibility in the continuum limit and still be consistent field theories. It is hard to imagine a more 
well-established quantum field theory than the two-dimensional 0(3) sigma model, so we would be 
reluctant to conclude that this model cannot be defined as the continuum limit of a lattice model. 
What seems more sensible is to regard the topological susceptibility as somewhat analogous to a 
mass parameter in a field theory with two or more independent mass scales. At least in more benign 
cases, sending one mass to infinity removes a sector of the theory but leaves the remainder of the 
theory well-defined. Such an interpretation of infinite topological susceptibility is also suggested by 
a Witten-Veneziano type argument relating the magnitude of xt to the splitting between the singlet 
and non-singlet mesons. By this argument, zz annihilation diagrams, allowed only in the singlet 
channel, will increase the mass 2 of the flavor singlet meson by an amount proportional to \t- The 
0(3) sigma model is an exactly solvable model, and it's particle spectrum is known to consist of 
only an isotriplet of mesons, with no isosinglet. From the quark model point of view, this suggests 
that the mass of the singlet has been shifted upward, if not by an infinite amount, then at least 
enough to make it unstable. By all indications, the mass shift of the flavor singlet meson is the only 
singular physical effect of the divergent xt, and the resulting theory is in fact well-defined. Direct 
spectroscopic studies of the singlet channel for various values of N should help to clarify this issue. 
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FIG. 12: Scaling behavior of the positive and negative contributions to the topological susceptibility in 
CP 1 and CP 3 
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The lowest dimension operator that couples to the singlet pseudoscalar channel is the topological 
charge q(x). As we have seen, the two-point correlator for q(x) is dominated by contact terms and 
has little if any contribution from propagating intermediate states. However, the coupling of q{x) 
to singlet zz states may be suppressed, so it is difficult to even put a lower bound on the singlet 
mass gap from the q{x) correlator data alone. The study of other operators which couple to the 
pseudoscalar singlet channel (e.g. smeared zz operators) is essential for a complete resolution of 
this issue. 

The fact that the usual scaling xt ~ ^ 2 is violated by the contribution of small instantons is 
not surprising, considering that in the limit of zero lattice spacing, a small instanton consists of a 
unit of topological charge within a vanishingly small radius, so that the topological charge density 
is becoming singular in the continuum limit. The topological charge is much more spread out in 
the case of domain wall excitations and should have a better chance of scaling properly. However, 
in both CP N_1 [i| and in QCD 1], the domain walls are found to be of order lattice spacing in 
thickness, so that, while they are less singular distributions than the small instantons, they are 
still not smooth on the lattice spacing scale and would seem to have the potential to cause some 
problems with the scaling of topological susceptibility. The naive argument that scaling occurs 
because everything becomes smooth and slowly varying at distances comparable to the lattice 
spacing is clearly not applicable here. The numerical result that, for CP 3 and higher, topological 
susceptibility scales like fj 2 and the physical xt approaches a finite value in the continuum limit 
is thus not immediatedly obvious and probably deserves a deeper explanation. We note that in 
calculating xt the integral under the positive peak from < x < 2 and the integral of the negative 
tail for x > 2 cancel with increasing accuracy as beta becomes large, as shown in Fig. |12(a)| and 



12(b) This figure illustrates that, for the larger N models, the positive and negative contributions 



to the correlator integral are separately diverging in the continuum limit, but their sum scales 
properly to a finite topological susceptibility. For the N > 3 models, this cancellation reflects the 
dominance of the 5"(t) term, whose coefficient scales to a constant, compared to the coefficient of 
the 8(t) term, which scales oc /i 2 . 

One of our main motivations for undertaking a study of topological charge structure in the 
(jpN-i moc i e ig was t try to explore the theoretical significance of the recently discovered 
membrane-like coherent topological charge structure in QCD Monte Carlo configurations 1]. As 
pointed out by Liischer Jjj, one can at least naively estimate the instanton melting point in 4- 
dimensional SU (N) gauge theory. In Section III we showed that with our lattice formulation of 
the model, the action of a small instanton was approximately e = ^ x 6.8 . . . The critical value 
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e = 4-7T gives an instanton melting point of N cr a = 3.7. It may be possible to increase the small 
instanton action and thus lower N cr n by improvement of the lattice CP N_1 action, but the lattice 
e is bounded above by the continuum instanton action e < 2ttN, so that any lattice action for the 
CP N ~ l model will give 

N crit > 2 (71) 

(Thus, in particular, no amount of improvement will completely eliminate small instantons from 
the CP 1 model.). The lowest possible value of N cr u = 2 can be estimated within the renormalized 
continuum theory alone, without reference to a lattice formulation. It corresponds to the "tipping 
point" in the integral over instanton radius in a semiclassical instanton calculation. This is where 
the integral changes from being divergent at the small instanton end to being divergent at the 
large instanton end. (For example, the integral over instanton radius in the CP 1 model goes like 
J dp/ p.) The analogous integral in 4-dimensional SU{N) gauge theory (using the one-loop beta 
function) behaves like 

dPfLNfi (72) 



which has it's tipping point at 



P 



N crit = ^ ( 73 ) 



By this estimate, the value N = 3 for real QCD is comfortably on the large- A^ side of the critical 
value, but (|73*|) is only a lower bound on the value of N cr a for any particular lattice action, so it 
is not clear how seriously to take this estimate. The Monte Carlo evidence jjj that membranes, 
not instantons, dominate the topological charge distribution in QCD adds support for the notion 
that N = 3 QCD is in the large-A" regime, above the instanton melting point. As is well-known, 
the integral over instanton radius in QCD is strongly weighted toward large instantons and must 
be cut off in an ad hoc way at about 0.3/m. From the present point of view this is an indication 



12], 



that instantons are unstable toward melting or unwinding. As Liischer originally argued 
by periodicity in 9, the vacuum inside a Wilson bag of charge k = 1 has the same energy as 
the 8 = vacuum indicating that the confining force between bag walls is completely screened 
by the appearance of an anti-bag. The analogous process in the CP N ~ l model is just ordinary 
string breaking and the vanishing of the area term for integer-charged Wilson loops due to charge 
screening. With the confining force screened, the bag-antibag dipole layer is free to expand and 
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produce extended membranes with layers of topological charge of opposite sign juxtaposed, just 
the type of structure that is seen in the Monte Carlo configurations. 
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